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Abstract. XMM-Newton archival observations of the Be/X-ray binary candidate RX J0103.6—7201 revealed pulsations with 
a period of ~1323 s. This makes RX J0103.6—7201 the X-ray pulsar with the longest period known in the Small Magellanic 
Cloud (SMC). More than 150 X-ray observations of RX J0103.6—7201 by ROSAT, Chandra and XMM-Newton show flux 
variations by a factor of 50 on time scales of days to years. Using the accurate positions obtained from ACIS-I images the 
optical counterpart is identified with a V = 14.6 mag emission line star. EPIC spectra of RX J0103.6—7201 above 1 keV are 
consistent with an absorbed power-law with column density between (6—9)x 10^*^ cm“^, except during one observation when 
an extraordinary high value of 1.1 x 10^® cm~^ was measured which strongly attenuated the power-law emission below 3 keV. 
A soft excess between 0.5 and 1.0 keV is evident in the spectra which becomes best visible in the highly absorbed spectrum. 
The soft component can be reproduced by a thermal plasma emission model with its luminosity strongly correlated with the 
total intrinsic source luminosity. Including results from three other SMC Be/X-ray binaries extends the linear correlation over 
three orders of magnitude in source intensity, strongly suggesting that the same mechanism is responsible for the generation of 
the soft emission. 

Key words, galaxies: individual: Small Magellanic Cloud - stars: neutron - X-rays: binaries - X-rays: galaxies - individual: 
RXJ0103.6-7201 


1. Introduction 

High mass X-ray binaries (HMXBs) are binary systems con¬ 
sisting of a compact stellar object accreting matter from a high 
mass early type star. In most cases the detection of pulsations 
in the X-ray flux with periods in the range of 0.07 s to 1500 s 
reveals a spinning neutron star orbiting an OB supergiant or 
a Be star. The large population of HMXBs makes the SMC 
unique among the galaxies of the Local Group. First, the abso¬ 
lute number of known HMXBs is similar to that in the Milky 
Way, despite the comparatively small mass of the SMC and sec¬ 
ond, while the ratio of supergiant to Be binaries in the Milky 
Way is 1/3.5 as found from a detailed literature search, the vast 
majority of HMXBs in the SMC has a Be star counterpart. 

Currently 48 X-ray pulsars are known in the SMC (Coe 
et al. 2005), from which 45 are most likely in HMXB systems 
as suggested by their X-ray properties (Haberl & Pietsch 2004). 
Optical identifications confirmed 33 Be counterparts while 
only one supergiant system (SMCX-1) is known. Additional 
candidates for Be/X-ray binaries were found by Haberl & 
Sasaki (2000) from correlations of ROSAT catalogues of SMC 
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sources with lists of Ha emission line objects (Meyssonnier & 
Azzopardi 1993, hereafter MA93) and by Sasaki et al. (2003). 


One of the Be/X-ray binary candidates is 
RXJ0103.6—7201, its ROSAT position (with an uncer¬ 
tainty of 4") is consistent with the emission line object 
[MA93] 1393 (Haberl & Sasaki 2000). This hard ROSAT 
source was detected in PSPC (Haberl et al. 2000) and HRl 
(Sasaki et al. 2000) pointed observations of the SMC. From 
three ROSAT observations a factor of three flux variability 
was seen on time scales of years. A first analysis of early 
XMM-Newton data was performed by Sasaki et al. (2003). 
RX J0103.6—7201 is located close to the supernova remnant 
IE 0102—72.3 in the North-East of the SMC which is fre¬ 
quently observed as calibration source for the instruments 
on-board XMM-Newton and Chandra. About 150 Chandra 
observations were performed with RX J0103.6—7201 located 
on various ACIS CCD chips which are readily available 
from the Chandra archive. In addition twelve XMM-Newton 
observations with the EPIC instruments in various readout 
modes makes RX J0103.6—7201 the most frequently observed 
extragalactic X-ray binary. We analyzed the temporal and 
spectral properties of RX J0103.6—7201 and present the 
results in this paper. 
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Fig. 1. Periodograms of the Zf test from two XMM-Newton observations of RX J0103.6—7201 (top). Corresponding pulse profiles folded with 
the most likely period for different energy bands, together with hardness ratios (HR) are shown in the bottom panels. For sinusoidal variations 
pulsed fractions in the 0.3—8.0 keV band are 15.7±1.7% (rev. 521) and 13.2±2.2% (rev. 616). 


2. Data analysis 

The Chandra ACIS calibration observations of IE0102—72.3 
are typically of 8 ks exposure (five observations have 20 ks) 
and were used to determine an accurate position (Sect. 2.4) 
and to investigate the flux history of RX J0103.6—7201 (Sect. 
2.3). The twelve XMM-Newton EPIC observations performed 
until 2004, Oct. 15 have longer exposures of 20—30 ks (only 
one exception with 13 ks). The EPIC-pn (Striider et al. 2001) 


and EPIC-MOS (Turner et al. 2001) data were, therefore, used 
for a spectral and temporal analysis, in particular to search for 
X-ray pulsations. During five of the EPIC observations (here 
as observation a single pointing within a satellite revolution 
is defined) two different CCD readout modes were used for 
the pn camera (switching between EE: full frame, LW; large 
window and SW; small window mode) while the MOS cam¬ 
eras were operated in LW mode during the whole observa¬ 
tion (only during revolution 433 and 447 part of the time was 
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spent in timing mode). Different optical light blocking filters 
(thin, medium and thick) were chosen in different revolutions 
(see also Table 1). During revolutions 616 and 803 the MOS 
LW mode held of view did not cover RX J0103.6—7201. Also 
the MOS timing mode observations cannot be used. The data 
were processed using the XMM-Newton analysis package SAS 
6.0.0 to produce the photon event hies and binned data prod¬ 
ucts such as images and spectra. Events for spectral and tem¬ 
poral analysis of RX J0103.6—7201 were extracted from circu¬ 
lar regions (radius 20") around the source position and from 
nearby source-free regions for background spectra. Photon ar¬ 
rival times were corrected to the solar system barycentre. For 
the purpose of searching for relatively long periods of sec¬ 
onds to tens of minutes as observed from Be/X-ray binaries the 
EPIC-pn data from one satellite revolution with different read¬ 
out modes were merged. MOS and pn data were combined for 
the common time intervals, i.e. MOS data at the beginning of 
an observation when the pn camera was still calculating the off¬ 
set map (and also when the readout mode was changed) were 
disregarded. 

2.1. X-ray pulsations 

We performed a timing analysis of the EPIC data of 
RX J0103.6—7201, searching for pulsations in the broad en¬ 
ergy band (0.3—8.0 keV) using the Rayleigh Zf technique in 
the range 1—3000 s. Strong peaks in the probability density 
function were found from the observations in revolutions 521 
(pn SWh-LW and MOSl/2 LW), 616 (pn SWh-LW), 721 (pn FF 
and MOS 1/2 LW) and 803 (pn LW) with Zj of 93.2, 77.9, 68.4 
and 29.9, respectively. The maximum Zf of 93.2 derived for 
revolution 521 corresponds to a period detection conhdence of 
1.0—1.3x10“^® (disregarding the detection in other observa¬ 
tions). Fig. 1 shows two example periodograms (top) together 
with the EPIC X-ray light curves folded with the most likely 
period value derived from the test. To investigate energy 
dependencies of the pulse prohle, light curves in the 0.3 —1.5 
keV (soft), 1.5—8.0 keV (hard) and the total 0.3—8.0 keV en¬ 
ergy bands were created. Hardness ratios derived by dividing 
the pulse prohles in the hard and soft bands show signihcant 
variations caused by weaker pulse modulations in the soft band 
compared to the hard band. This is particularly evident during 
the revolution 616 observation where no signihcant modula¬ 
tion was seen in the soft band. Period errors (Icr) were deter¬ 
mined by htting of a sine wave to the total light curves with 
80 s binning, taking into account long term trends visible in 
the light curves by adding polynomial components to the sine 
model light curve. We also tested the inclusion of a second har¬ 
monic component as might be suggested by the analysis shown 
in Fig. 1. This did not improve the results because the major 
uncertainties in the determination of the long pulse period arise 
from the superposed long term trends. The pulse period his¬ 
tory for the four observations is shown in Fig. 2. No signihcant 
change in the pulse period is visible consistent with a constant 
period of (1323.2 ± 5.8) s. However, the large errors are still 
compatible with large spin-up or spin-down with up to ± 13 s 
year“^. 
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Fig. 2. Pulse period history of RX J0103.6—7201 over 1.5 years. The 
pulse periods of 1322.1±5.9, 1327.0±6.0, 1319.3±7.5 and 1322±12 
were derived from the EPIC data from revolutions 521, 616, 721 and 
803, respectively. 


2.2. X-ray spectra 

Pulse phase averaged EPIC-pn and -MOS spectra were ex¬ 
tracted using pattern 0—4 (single-tdouble pixel) and pattern 
0—12 (single — quadruple) events, respectively, and binned to 
obtain at least 50 counts per bin. For observations split into 
parts with different readout modes (with different spectral re¬ 
sponse) separate spectra were obtained. Nine EPIC observa¬ 
tions yielded a sufficiently high number of counts for a spec¬ 
tral analysis. All spectra for each observation were ht simul¬ 
taneously (using XSPEC vll.3) with the same model only al¬ 
lowing a free normalization factor between the different instru¬ 
ments and readout modes (which cover different time intervals 
with varying source intensity). Errors were determined for 90% 
conhdence levels. A simple power-law model including photo¬ 
electric absorption by matter with solar abundance reproduces 
the data well above 1 keV while at lower energies a soft excess 
is apparent, in particular in the EPIC-pn spectra. This soft com¬ 
ponent is clearly visible during the XMM-Newton observation 
of revolution 888 when the power-law component was highly 
absorbed and did not contribute below 3 keV. The EPIC spec¬ 
tra of this observation were used to investigate the nature of the 
soft component. 

To ht the soft component a black-body model and a model 
for a hot plasma in collisional equilibrium (Mewe et al. 1985, 
MEKAL model in XSPEC) were compared. These were used 
in the past to reproduce low energy excesses seen in spectra of 
SMC Be/X-ray binaries (Sasaki et al. 2003; Haberl et al. 2003; 
Majid et al. 2004). The EPIC spectra of RXJ0103.6-7201 
from revolution 888 were ht with these models between 0.2 and 
2.0 keV where the power-law contribution can be neglected. 
This results in x^/dof values of 19.04/15 and 14.50/12 for the 
black-body and the thermal model, respectively. The thermal 
plasma model better reproduces the structures in the spectrum 
which are probably caused by emission lines around 650 eV 
and 850 eV. Therefore, the thermal model was used as addi¬ 
tional component to the power-law to ht the full EPIC spectra 
from all observations. 
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Fig. 3. EPIC spectra of RX J0103.6—7201 from six observations. 


For the revolution 888 spectra two individual column den¬ 
sity values were used for the two spectral components result¬ 
ing in very different best ht values for Nh of (5. 7^2 3 )^ 10^^ 
cm“^ and 107±28 x 10^^ cm“^ for the thermal and the power- 
law component, respectively. For all other observations the Nh 
values were not significantly different and therefore the same 
(but free to vary in the fit) value was used for both components. 
Also the temperature was fixed at the value derived from the 
revolution 888 observation (kT = 0.15±0.04 keV). In Table 1 
the results for the spectral fits are summarized; photon index 
7, absorption column density Nh, emission measure EM = f 
nnUe dV, together with fluxes and source luminosities and re¬ 
duced values per degree of freedom (dof). In Fig. 3 the EPIC 
spectra together with the best fit models are plotted. 


2.3. Energy dependence of pulsations 


The pulsed fraction observed from RXJ0103.6—7201 dur¬ 
ing revolution 616 is largely reduced (if present at all) in 
the 0.3—1.5 keV band compared to the 1.5—8.0 keV band. 
However, the observation from revolution 521 shows that there 
can be still a large contribution from the power-law component 
in the soft band. The observation from revolution 888 where 
the soft component was completely isolated is ideally suited to 
investigate the pulse profiles of both spectral components sep¬ 
arately. This is shown in Eig. 4 and although the statistics is 
relatively poor, pulsations can be detected above 1.5 keV while 
below they are strongly reduced or may not be present at all. 
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Table 1. XMM-Newton EPIC spectral fit results. 


Instr. mode 

Rate’-'^^ 

Nh 

7 

EM 

Flux^^^ 

T (2) 

T i (3) 

X?/dof 

filter 

cts 

lO^^cm-^ 

10®® cm-® 

erg cm'® 

erg 

erg 

2000 Apr. 16/17 (Revolution 065) 







PN FF thin 

0.177 

6.5±L5 

0.88±0.07 

6.93 

1.61x10“^^ 

6.87x10®® 

2.00x10®® 


Ml LWthin 

0.052 



7.18 

1.66x10"^^ 

7.13x10®® 

2.08x10®® 


M2 LW thin 

0.051 



6.90 

1.60x10-^^ 

6.84x10®® 

2.00x10®® 


PN FF medium 

0.156 



6.10 

1.41x10"^^ 

6.04x10®® 

1.76x10®® 


Ml LW medium 

0.048 



6.56 

1.52x10-^^ 

6.50x10®® 

1.90x10®® 


M2 LW medium 

0.049 



6.63 

1.54x10“^^ 

6.57x10®® 

1.92x10®® 

1.29/209 

2001 Apr. 14/15 (247) 








PN SW thin 

0.413 

6.2±0.4 

0.84±0.03 

22.1 

3.59x10“^^ 

1.54x10®® 

5.61x10®® 


PN FF thin 

0.376 



20.0 

3.25x10"^^ 

1.39x10®® 

5.08x10®® 


Ml LWthin 

0.121 



21.9 

3.55x10-^^ 

1.52x10®® 

5.55x10®® 


M2 LW thin 

0.126 



21.9 

3.57x10"^^ 

1.53x10®® 

5.58x10®® 

1.24/347 

2001 Dec. 25/26 (375) 








PN SW thin 

0.189 

7.0±L6 

0.90±0.08 

6.03 

1.77x10“^^ 

7.56x10®® 

1.94x10®® 


Ml LWthin 

0.064 



6.86 

2.01x10-^^ 

8.61x10®® 

2.20x10®® 


M2 LW thin 

0.063 



6.73 

1.97x10“^^ 

8.44x10®® 

2.16x10®® 

1.19/177 

2002 Oct. 13 (521) 
PN SW thin 

0.244 

6.8±L0 

0.85±0.05 

12.1 

2.35x10"^^ 

1.01x10®® 

3.27x10®® 


PN LW thin 

0.289 



14.7 

2.88x10"^^ 

1.23x10®® 

4.00x10®® 


Ml LWthin 

0.096 



15.4 

3.00x10"^^ 

1.28x10®® 

4.16x10®® 


M2 LW thin 

0.094 



15.0 

2.91x10“^^ 

1.24x10®® 

4.04x10®® 

1.08/185 

2002 Dec. 14 (552) 
PN SW thin 

0.244 

6.5±L3 

0.91±0.09 

4.83 

5.50x10“^® 

2.36x10®® 

1.12x10®® 


PN LW thin 

0.289 



5.74 

6.56x10“^® 

2.81x10®® 

1.33x10®® 


Ml LWthin 

0.096 



5.63 

6.41x10-^® 

2.75x10®® 

1.30x10®® 


M2 LW thin 

0.094 



5.83 

6.64x10"^® 

2.84x10®® 

1.35x10®® 

0.82/63 

2003 Apr. 20/21 (616) 








PN SW medium 

0.200 

8.4±L2 

0.96±0.06 

30.9 

1.88x10"^^ 

8.04x10®® 

6.40x10®® 


PN LW medium 

0.166 



26.9 

1.63x10-^^ 

6.98x10®® 

5.55x10®® 

0.99/111 

2003 Nov. 16 (721) 
PN FF thick 

0.143 

6.2±L1 

0.91±0.05 

8.53 

1.32x10-^^ 

5.65x10®® 

2.14x10®® 


Ml LWthin 

0.050 



9.55 

1.48x10“^^ 

6.32x10®® 

2.40x10®® 


M2 LW thin^^) 

0.036 



7.66 

1.18x10“^^ 

5.07x10®® 

1.92x10®® 

1.48/125 

2004 Apr. 28 (803) 
PN LW thick 

0.030 

6.6±L7 

L12±0.15 

3.28 

2.82x10“^® 

1.21x10®® 

7.21x10®® 

0.72/25 

2004 Oct. 14 (888)^®) 








PN LW thick 

0.055 

107±28 

L03±0.34 

4.19 

7.31x10“^® 

3.13x10®® 

1.34x10®® 


Ml LW thick 

0.013 



4.40 

7.67x10“^® 

3.28x10®® 

1.40x10®® 


M2 LW thick^^^ 

0.008 



3.18 

5.56x10“^® 

2.38x10®® 

1.02x10®® 

1.48/45 


Net count rates and observed total fluxes in the 0.2—10.0 keV band. 

X-ray 0.2—10.0 keV luminosity (including absorption) for a distance of 60 kpc. 

Source intrinsic X-ray luminosity in the 0.2—10.0 keV band (corrected for absorption). 

Source on bad CCD column. 

Individual column density for the thermal plasma component Nh = (5.7j(2'3)x 10^^ cm~^. The temperature of the thermal plasma model 
component was determined from this observation to kT = 0.15±0.04 keV and fixed in the fits to the other observations. 


2.4. Long-term intensity variations 

To look at intensity variations of RX J0103.6—7201 on time 
scales of days to years we converted ROSAT and Chandra 
count rates obtained from source detection analyses into 
0.2—10.0 keV fluxes. We assumed a power-law plus MEKAL 
spectrum with photon index 0.9, column density of 6.5x10^^ 


cm“^ and temperature 0.15 keV, typical values as derived from 
the EPIC spectra. Similarly, for EPIC observations with insuffi¬ 
cient statistical quality (short exposure and/or low source flux) 
count-rates were used. Eor the estimates of ROSAT fluxes rel¬ 
atively large systematic errors may be introduced due to the 
strong influence of the assumed column density on the soft 
0.1—2.4 keV band count rate to luminosity conversion (an ad- 
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Fig. 4. EPIC-pn light-curve from revolution 888 in soft and hard en¬ 
ergy bands folded with the pulse period. 

ditional Nh of 1.0x10^^ cm~^ changes the conversion factor 
by 12%). The measured fluxes as function of time over nearly 
13 years are plotted in Fig. 5. There seems to be a base flux 
level of ~7.6xl0“^^ erg cm“^ s“^ below which the source 
was never seen. The maximum flux of 3.7xl0~^^ erg cm“^ 
s“^ is a factor of 50 higher and was reached during the XMM- 
Newton observation of revolution 247 (second part with EPIC- 
pn in FF mode). A similar outburst level was seen also dur¬ 
ing the third ROSAT observation, four Chandra and one other 
XMM-Newton observation. Chandra observations were often 
performed in blocks spread over ~0.7 days during which the 
source stayed at a similar intensity level, indicating that high 
and low states last longer than that. The source was generally 
brighter (both, low and high intensity levels decreased) dur¬ 
ing the first half as compared to the second half of the last five 
years (Fig. 5) which may be related to long-term changes in the 
properties of the circum-stellar matter around the mass donor 
star. 

2.5. Source identification 

X-ray source positions were obtained from the analysis of 
all observations where the source was detected in the ACIS-I 
CCDs. The average position from 71 observations with a stan¬ 
dard deviation of I'.'l is only O'.'07 away from the 2MASS ob¬ 
ject 01033752-7201330 at RA = 0l‘'03“37!524 and Dec = 
-72°01'33'.'03 (J2000.0) for which an error of O'.'ll is given 
(Cutri et al. 2003). For comparison the average EPIC position 
is within O'.'S of the 2MASS position. This star is covered by 
the SMC UBVR CCD survey of Massey (2002) and also listed 
in the UBVI photometry catalogue of Zaritsky et al. (2002). 
Optical positions and available magnitudes are summarized in 
Table 2. The star can be identified with the emission line object 
([MA93] 1393, Meyssonnier & Azzopardi 1993). The presence 
of a ~14.6 Vmag star inside the error circle of the Chandra 
position clearly establishes it as the optical counterpart of the 
X-ray source. 

The inferred effective temperature and bolometric luminos¬ 
ity of the optical counterpart given by Massey (2002) suggest 



4.9x10* 5x10* 5.1x10* 5.3x10* 5.3x10* 

MJD (days) 


Fig. 5. X-ray flux variations from RX J0103.6—7201 since the first 
ROSAT (triangles) observation in Nov./Dec. 1991. Circles mark 
XMM-Newton observations with fluxes directly derived from spec¬ 
tral fits to the EPIC-pn spectra (which yield two data points when two 
read-out modes were used. All other fluxes were converted from count 
rates assuming a power-law -l- MEKAL model with typical parameters 
as derived from the EPIC spectra (see text). Crosses indicate Chandra 
observations and squares XMM-Newton observation with the source 
at low intensity. 

a spectral type of 05 V. Assuming this spectral type the mea¬ 
sured B—V index imply E(B—V) of 0.32 and Nh = 1.9x10^^ 

_2 

cm 

3. Discussion 

3.1. RX J0103.6- 7201, a Be/X-ray binary puisar 

The discovery of a 1323 s modulation in the X-ray flux from 
RX J0103.6—7201 increases the number of pulsars in the SMC 
to 49 from which up to 46 are most likely neutron stars in 
HMXB systems. The pulse period of 1323 s is the longest 
known from SMC pulsars. It considerably exceeds that of 
AX J0049.4-7323 (= RX J0049.7-7323) (756 s; Yokogawa 
et al. 2003) and is similar to that of RX J0146.9-f6121 (Haberl 
et al. 1998) in the Milky Way. Optical spectra of 33 counter¬ 
parts from HMXBs in the SMC identify the majority as Be-star 
binaries (Coe et al. 2005), while only one (SMC X-1) is a super¬ 
giant system. Brightness and colours of the optical counterpart 
and the presence of Ha emission (MA93) strongly suggest that 
also RX J0103.6—7201 is a Be/X-ray binary system. 

The X-ray spectrum of RX J0103.6—7201 is typical for 
Be/X-ray binaries in the SMC. The power-law photon index 
of xoO.9 is well within the relatively narrow range observed 
from these systems (Haberl & Pietsch 2004). The absorbing 
column density was found to be similar during most obser¬ 
vations yielding (6—9)xl0^^ cm“^, except during the XMM- 
Newton observation in revolution 888 when Nh exceeded 10^^ 
cm“^. Such high values are usually seen only from supergiant 
HMXBs when the X-ray source is viewed through the dense 
innermost parts of the stellar wind of the supergiant before or 
after eclipse (e.g.4U1700—37 and Vela X-1; Haberl et al. 1989; 
Sako et al. 1999). This suggests that we may view the system of 
RX J0103.6—7201 at high inclination and during 2004, Oct. 14 
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Table 2. Optical identifications. 


Colours (catalogue) 

R.A. and Dec. (2000.0) 

Vmag 

B-V 

U-B V-R V-I 

UBVR (Massey 2002) 

UBVI (Zaritsky et al. 2002) 

0l‘'03™37!46 -72°0T33'.'2 
01®03™37;52 -72°0T33'.'0 

14.54 

14.65 

-0.06 

-O.ll 

-0.99 -0.08 

-1.08 - -0.04 


(revolution 888) the line of sight to the X-ray source crossed the 
wind region very close to the Be star. Since Be/X-ray binaries 
are much wider systems in comparison to supergiant systems 
such cases should be very rare. Indeed, RX J0103.6—7201 is 
to our knowledge the only Be/X-ray binary which has shown a 
spectrum with such extreme absorption. Alternatively absorp¬ 
tion by dense material close to the neutron star cannot be ruled 
out, but may be expected to occur more often. 

3.2. Soft emission in the X-ray spectra of Be/X-ray 
binaries 

A soft spectral component in the X-ray spectrum of 
RX J0103.6—7201 was first reported by Sasaki et al. (2003) 
who analyzed the XMM-Newton data from revolution 247 
and fit the EPIC-pn spectrum with thermal plasma and power- 
law components. This soft emission became clearly visible in 
the revolution 888 observation when both spectral components 
were separated into different energy bands. From Table 1 it is 
evident that the intensity of the soft component, given as emis¬ 
sion measure (which directly reflects the intrinsic luminosity 
of the soft component) in the range of (3—30)x 10^® cm ® 
strongly correlates with the intrinsic source luminosity. This is 
demonstrated in the top panel of Fig. 6 where the luminosity of 
the soft component is plotted versus the total source luminosity. 

Similar soft emission was also observed from the 
pulsars EXO 053109-6609.2 with XMM-Newton (Haberl 
et al. 2003) and RX J0059.2—7138 in a simultaneous 
ROSAT/ASCA observation (Kohno et al. 2000) with emis¬ 
sion measures of ~4xl0®^ cm“® and ~10®^ cm“®, respec¬ 
tively. RX J0059.2—7138 was observed at a luminosity of 
2.5x10®® erg s“®, i.e. both luminosity and emission mea¬ 
sure were a factor of 40 higher compared to the maxi¬ 
mum values observed from RX J0103.6—7201. On the other 
hand, for EXO 053109—6609.2 an intrinsic source luminos¬ 
ity of 4.6x10®^ erg s“® was reported which yields a ratio 
EM[cm“®]/Lx[erg s“®] of 8.7x10^® a factor of 20 higher 
than for RX J0103.6—7201. Sasaki et al. (2003) report two 
more SMC pulsars which exhibit a soft emission component. 
We re-analyzed the EPIC-pn spectra of RXJOIOI.3—7211 
and AXJ0103—722 in a consistent way as was done for 
RX J0103.6—7201. Keeping the temperature of the thermal 
plasma component fixed at 0.15 keV and using a common Nh 
for soft and hard (power-law) component we derive EM and in¬ 
trinsic source luminosity (0.2—10.0 keV). The results, together 
with the values found in the literature, are included in the bot¬ 
tom panel of Fig. 6. They are consistent with the linear rela¬ 
tion between EM (or luminosity of the soft component) and to¬ 
tal source luminosity derived for RX J0103.6—7201 with EM 


[cm“®] = 4x 10^^ xLx [erg s“®] (plotted as straight line). Here 
we use EM instead of luminosity of the soft component to be 
able to compare with published data. 

Hickox et al. (2004) explored four possible mechanisms for 
the origin of the soft excess in X-ray pulsars: 1) re-processing 
of hard X-rays by optically thick, dense material near the neu¬ 
tron star; 2) re-processing by optically thin material; 3) thermal 
emission from collisionally ionized plasma and 4) soft emis¬ 
sion from the accretion column and the neutron star surface. 
They conclude that for luminous sources (Lx > 10®® erg s ®) 
only case 1) can explain the strong soft component while for 
luminosities below 10®® erg s“® all three other cases can con¬ 
tribute. 

The tight correlation between the luminosity of the soft 
spectral component and the total source luminosity we found 
for RX J0103.6—7201 over a factor of 10 in luminosity and 
moreover extending to three orders of magnitude when in¬ 
cluding results from three other SMC Be/X-ray binary pulsars 
strongly suggests that there is a single mechanism at work over 
a wide range of source intensity. The correlation as such makes 
emission from a collisionally ionized plasma unlikely and sup¬ 
ports the models involving re-processing of the X-ray emission 
from the neutron star. Therefore, the MEKAL model is actu¬ 
ally not the proper model to derive parameters of the emitting 
plasma like temperatures or elemental abundances. However, 
it best fits the soft part of the EPIC spectra among all avail¬ 
able simple XSPEC models and allows to accurately determine 
fluxes. 

The pulsed fraction observed from RX J0103.6—7201 is 
largely reduced (if present at all) at energies below 1.5 keV. 
Similarly, no pulsations at energies below 0.4 keV were seen 
from the 2.76 s pulsar RX J0059.2-7138 (Kohno et al. 2000) 
which indicates a size for a re-processing region with r > 
cxPspin ~ 10®® cm. To smear out longer pulse periods the 
size must be correspondingly larger and for RX J0103.6—7201 
would require a size > 4x 10®® cm. An independent very rough 
estimate for the size of the re-processing region can be derived 
from the emission measure f nnUe dV ~ 47r/3 r®xn^ and the 
absorbing column density f nn dr rxNn- For EM in the 
range of (3—30) x 10®® cm“® and a typical local Nh of 6x10^® 
cm“^, r is in the range 2x10®^ cm to 2x10®® cm (increasing 
with higher EM and hence source luminosity). If the large ab¬ 
sorption observed from revolution 888 for the power-law com¬ 
ponent is caused by dense stellar matter close to the Be star, 
than these latter size estimates must be regarded as upper limit. 
A lower limit of 1 x 10®^ cm can then be derived using the high 
Nh value. These simple estimates argue for reprocessing of the 
hard X-rays in a large part of the optically thin stellar wind. The 
correlation of soft flux with source intensity for four Be/X-ray 
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Fig. 6. Top: Luminosity of the soft spectral component (modeled as 
MEKAL component) in the EPIC spectra of RX J0103.6—7201 vs. 
total source luminosity. Both luminosities are source intrinsic, i.e. cor¬ 
rected for photo-electric absorption. Bottom: Emission measure of the 
thermal plasma component as function of the total source intrinsic lu¬ 
minosity including results from other Be/X-ray binaries. 
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binaries in the SMC suggests that this mechanism is the dom¬ 
inant process to create the soft spectral component at least in 
HMXBs with a Be star as mass donor. 

It remains to be explained why EXO 053109—6609.2 ex¬ 
hibits a much stronger soft component compared to the SMC 
Be/X-ray pulsars. It is unlikely that the higher metallicity of 
the EMC can cause a factor of 20 increase in reprocessing 
efficiency. EXO 053109—6609.2 has shown a more complex 
spectrum during the XMM-Newton observation reported by 
Haberl et al. (2003) which may be interpreted as partial cov¬ 
ering of the X-ray source. Hence, a large fraction of the hard 
emission might be completly blocked and the total source lu¬ 
minosity underestimated. Other Be/X-ray binaries in the EMC 
with a soft emission component need to be observed to verify if 
EXO053109—6609.2is an exceptional case (e.g. A0538—66, 
Mavromatakis & Haberl 1993). 
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